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Abstract. A growing interface subject to noise is described by the Kardar- 
Parisi-Zhang equation or, equivalently, the noisy Burgers equation. In one 
dimension this equation is analyzed by means of a weak noise canonical 
phase space approach applied to the associated Fokker-Planck equation. The 
growth morphology is characterized by a gas of nonlinear soliton modes with 
superimposed linear diffusive modes. We also discuss the ensuing scaling 
properties. 
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1. Introduction 

Macroscopic phenomena far from equilibrium abound and include phenomena such as 
turbulence in fluids, interface and growth problems, chemical reactions, processes in 
glasses and amorphous systems, biological processes, and even aspects of economical 
and sociological structures. 

In recent years much of the focus of modern statistical physics and soft condensed 
matter has shifted towards such systems. Drawing on the case of static and dynamic 
critical phenomena in and close to equilibrium, where scaling, critical exponents, and 
the concept of universality have so successfully served to organize our understanding 
and to provide a variety of calculational tools, a similar approach has been advanced 
towards the much larger class of nonequilibrium phenomena with the purpose of 
elucidating scaling properties and more generally the morphology or pattern formation 
in a driven state. 

In this context the noisy Burgers equation or the equivalent Kardar-Parisi-Zhang 
(KPZ) equation, describing the nonequilibrium growth of a noise-driven interface, 
provide simple continuum models of an open driven nonlinear system exhibiting scaling 
and pattern formation. 

In one dimension the KPZ equation has the form (in a co-moving frame) ^ 

^=i^V'h+^{Vhr+r,. (1) 

Here h is the height field, ly a damping or viscosity characterizing the linear diffusive 
term, A a coupling strength for the nonlinear mode coupling or growth term, and 
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Figure 1. We depict the growth morphology of a growing interface in a system 
of size L. The saturation width in the noise-driven stationary state is denoted by 

UIQ. 




Figure 2. We depict the interface width w{t) as a function of t. In the transient 
regime for t <^ to ^ w grows according to t^/^. In the stationary regime 
attained for t ^ to the width w saturates to the value wo ^ ■ 

rj a Gaussian white noise, driving the system into a stationary state. The noise is 
correlated according to 

{ii{xt)T^m) = ^s{x)m , (2) 

and characterized by the noise strength A. In figure ^ we have depicted a reahzation 
of a growing interface. 

Assuming an initiaUy flat interface the width w{t, L) grows in time approaching 
the saturation width wq in the stationary regime. The dynamical scaling hypothesis 
Q then asserts that 

w{t,L)^L^F^{t/L') , (3) 

where C is the roughness exponent characterizing the morphology of the interface and 
z the dynamic exponent describing the dynamical correlations. The scaling function 
F^(u) together with C and z define the scaling universality class and has the limits 
F^(u) — > const, for u ^ oo and F„{u) — + u^/^ for u 0. In figure ^ we have shown 
the time- and size dependent width 'w{t, L) Similarly, applying the dynamical scaling 
hypothesis to the height correlations in the stationary regime we have 

{{h{xt) - {h{xt))){hm - {h{QO)))) = x^<Fi,{t/x') , (4) 

where the scaling function Fh obeys Fh(u) — > const, for u ^ oo and Fh(u) — > u^^/^ 
for u 0. 

In what follows it turns out that the local slope field u = V/i is the appropriate 
variable in which to discuss the growth morphology of an interface. The height field 
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Figure 3. We depict the morphology of a growing interface in terms of the local 
slope. The slope field fluctuates about zero. 



h = J dx u is then an integrated variable sampling the slope fluctuations. In figure ^ 
we have depicted a realization of the slope field The dynamical scaling hypothesis for 
the slope field then reads 

{u{xt)u{00)) = x^'^-^F^{t/x'') , (5) 
with limits F^{u) const, for m ^ oo and F^iu) ■(i(2C-2)/z j-qj. ^ ^ q_ 

The stochastic dynamics of the slope field is then according to Eq. (|l]) governed 
by the noisy Burgers equation |^ 

du ^ 

— = jyV^u + AmVu + V77 . (6) 

This equation in the noiseless case for 77 = was originally proposed by Burgers 
in order to model turbulence in fluids; note the similarity with the Navier-Stokes 
equation for A = — 1. 

The substantial conceptual problems encountered in nonequilibrium physics are 
in many ways embodied in the KPZ-Burgers equations and describing the 
self-affine growth of an interface subject to annealed noise arising from fluctuations 
in the drive or in the environment, and these equations in one and higher dimensions 
have been the subject of intense scrutiny in recent years owing to their paradigmatic 
significance within a field theory of nonequilibrium systems ^, ||, |lO[ ^ ^ |l^, 

Interestingly, the Burgers-KPZ equations are also encountered in a variety of other 
problems such as randomly stirred fluids, dissipative transport in a driven lattice gas, 
the propagation of flame fronts, the sine-Gordon equation, and magnetic flux lines 
in superconductors. Furthermore, by means of the Cole-Hopf transformation the 
Burgers-KPZ equations are also related to the problem of a directed polymer or a 
quantum particle in a random medium and thus to the theory of spin glasses. 

In a series of papers the one dimensional case defined by (^) has been analyzed 
in an attempt to uncover the physical mechanisms underlying the pattern formation 
and scaling behavior. Emphasizing that the noise strength A constitutes the relevant 
nonperturbative parameter driving the system into a statistically stationary state, the 
method was initially based on a weak noise saddle point approximation to the Martin- 
Siggia-Rose functional formulation of the noisy Burgers equation [0 |l^, . This 

work was a continuation of earlier work based on the mapping of a solid-on-solid model 
onto a continuum spin model [ pT| . More recently the functional approach has been 
superseded by a canonical phase space method deriving from the canonical structure of 
the Fokker-Planck equation associated with the Burgers equation |3[ |^ In 
the present context we attempt to give a brief account of this approach with emphasis 
on the physical interpretation. 
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2. The noisy Burgers equation 

The noisy Burgers equation has the form of a conserved nonUnear Langevin 
equation, du/dt = — Vj, with fluctuating current j = —v\7u — (A/2)u^ + 77. The 
hydrodynamical origin of the Burgers equation, as reflected by the presence of the 
mode couphng or convective term AmVu, imphes that the Burgers equation is invariant 
under the Gahlean transformation 

X ^ X — XuQt and u u + uq , (7) 

involving a shift of the slope field. Since the nonlinear coupling strength A enters as a 
structural constant in the symmetry group it is invariant under scaling and a simple 
Kadanoff-type block renomalization group argument in both space and time implies 
the dynamical scaling law 

C + z = 2 , (8) 

providing a relationship between the roughness exponent and the dynamic exponent. 

In the linear case for A = the Burgers equation takes the Edwards- Wilkinson 
(EW) form (2§1 

du/dt ^ i^V^u + Vri , (9) 

i.e., a linear diffusion equation driven by conserved noise. This equation is easily 
decomposed and analyzed in terms of wavenumber modes Uk = J dx u{x) exp {—ikx), 
u^. = u-k- In the noiseless case for 77 = the field Uk, governed by an overdamped 
oscillator equation, decays according to Uk oc exp(— i/fc^t); Tk — \/vk^ setting the 
spectrum of relaxation times. For r/ ^ the slope field Uk is driven into a noisy 
stationary state. Defining Ukui = J dt Uk e'xp{iujt) the stationary correlations are given 

by 

We note that ( |l0| ) in addition to the pole oj — iiyk^, corresponding to the decaying 
mode, also has a pole at uj ~ —ivk"^, characterizing a growing mode. The noise in 
driving the equation thus excites a growing diffusive mode and we have generally for 
a particular realization 

u{xt) = {Ae-"^^'' + 56+"'='*) e*'^^ . (11) 

In the stationary state the correlation are time reversal invariant, requiring both 
growing and decaying modes; this also follows from the Fokker-Planck analysis in 
section 3. From (|lO|) we also obtain the static correlations (m(x)u(O)) = /S./2v5{x)^ 
showing that u{x) is spatially uncorrelated. Moreover, the stationary distribution has 
the form 



Pst(") OC exp 



dx u{xy 



(12) 



V 

'a 

Comparing (|l^) with the scaling form (^ we also infer the scaling exponents 
( = 1/2 and z — 2, characteristic of diffusion and defining the EW universality 
class. Also, noting that the diffusive term in (|^) can be derived from a free energy 
F = (1/2) J dxu^, it follows that the EW equation describes the fluctuations in an 
equilibrium system with temperature A/2i/, i.e., Pg^^ — exp [—{2h'/A)F]. 

Summarizing, in the linear case the elementary excitations are decaying and 
growing nonpropagating diffusive modes. Since from (|l|) d{h)/dt = vV'^{h) the 
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Figure 4. We show the slope and height fields for a constant shift no. The 
sloped height profile then moves with velocity Xuq. 




Figure 5. We show a step profile in the slope field with amplitudes ni and 
U2. The step moves with mean velocity v = {vi + V2)/2, where vi = Xui and 
V2 = \u2 are the velocities of the plateaus. The mode correspond to the slope 
dependent velocity of the height field subject to growth. 



average height (h) decays to zero (with respect to the co-moving frame) and the 
diffusive modes do not represent nonequiUbrium growth but in fact characterizes 
equihbrium fluctuations at the temperature A/2i'. FinaUy, the exponents C = 1/2 
and z = 2 for the EW universahty class reflect the uncorrelated slope field, or the 
random walk of the integrated slope, i.e., the height h, and the diffusional character 
of the dynamics, respectively. 

In the nonlinear case for A 7^ the Galilean invariance (^ becomes operational 
and we can draw some simple conclusions concerning a growing interface. First, 
imposing a constant shift uq (|^) implies that the tilted height field propagates with 
velocity Xuo as shown in figure ^ Next, considering a step in the slope field with 
amplitudes ui and U2 propagating with velocities vi = Aui and V2 = Aw2, it follows 
that the step itself propagates with mean velocity v = {vi+V2)/2. For the height field 
this corresponds to the part with the largest slope moving with the largest velocity 
in accordance with the slope dependent growth velocity in the KPZ equation (|^), i.e., 
d{h)/dt = (A/2)((V/i)^) > 0. That kind of localized modes giving rise to growth 
are indeed supported by the Burgers equation which possesses a spectrum of soliton 
modes. The configuration is depicted in figure ^ 

In the noiseless case for 77 = the Burgers equation support a right hand single 
soliton solution, in the static case of the kink-like form 

uix) = uotanh — -{x — xq) , (13) 
2v 

localized at xq, with amplitude uq, and width 2v/\uq. Boosting the soliton to a finite 
velocity this mode corresponds to the configuration depicted in figure |^. Denoting the 
right and left hand boundaries by uj^ and u_, respectively, we also infer the soliton 
condition 



u+ + U- — —2v/\ 



(14) 
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The relaxing growth morphology of the noiseless Burgers equation corresponding to 
the deterministic transient growth of the KPZ equation can be described by a gas of 
propagating right hand solitons, connected by constant-slope ramp solutions, and with 
a spectrum of superimposed decaying linear modes. In the height field this morphology 
corresponds to downward cusps connected by parabolic segments with superimposed 
linear modes |, |2|,|2|. 

In the noisy case the interface is driven into a stationary state, and anticipating 
the analysis in section 3, it turns out that the noise excites a left hand soliton of the 
shape 

u{x) = —uotetnh -^{x — xq) ■ (15) 

This mode is a solution of the growing noiseless Burgers equation for v —v. This 
doubling of modes is equivalent to the linear case where the mode has the form (pi]). 
The stationary growth morphology can thus be described by a gas of right hand and 
left hand solitons matched by the soliton condition ( p^ with superposed linear modes. 
The Fokker-Planck equation for the Burgers equation has the form 

Af^i^P, (16) 

where P is the probability distribution and H a Hamiltonian (Liouvillean) driving the 
equation of the form 

dx—iv^'^u + AuVu) 
bu 

+ -A^ / dxdx' VV'6{x - x')—— . (17) 
2 J ouou 

Whereas (16 will be discussed in more detail in section 3 it follows easily that the 
stationary solution of (|6|) has the form in (|l2) |2^, independent of A, and we infer 
as in the linear case C, = 1/2. The scaling law (8|) then implies the dynamic exponent 
z = 3/2. 

Summarizing, in the Burgers case the elementary excitations are propagating right 
hand and left hand solitons with superimposed linear modes. The soliton propagation 
corresponds to nonequilibrium growth. Finally, the scaling exponents C, = 1/2 and 
z = 3/2 defining the Burgers/KPZ universality class correspond to the uncorrelated 
slope field or random walk of h and soliton propagation, respectively. 



3. The weak noise approach 

Apart from numerical modeling and analysis of other models falling in the same 
universality class, the standard analytical approaches to the noisy Burgers equation 
are i) the dynamical renormalization group method (DRG) and ii) the mode coupling 
approach (MC). The DRG method accesses the scaling regime and provide an epsilon 
expansion about d — 2. Above d = 2 the system exhibits a kinetic phase transition 
from a smooth phase with EW exponents to a rough phase controlled by a strong 
coupling fixed point with largely unknown scaling exponents. In d = 1 the scaling 
is controlled by a strong coupling fixed point and the DRG yields (fortuitously) 
the known scaling exponents. Unlike the success of the DRG in dynamical critical 
phenomena, the results obtained for the noisy Burgers equation are limited despite 
a substantial theoretical effort. The MC method by neglecting (unrenormalized) 
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vertex corrections provides closed equations yielding scaling functions. However, the 
somewhat ad hoc nature of the MC approach makes it difficult to make contact with 
more systematic approaches. 

The functional or the equivalent phase space approach valid in the weak noise 
limit, A ^ 0, replaces the stochastic Langevin-type Burgers equation by coupled 
deterministic diffusion-advection type mean field equations, 
du 

— = j/V^M - V^p + AmVm , (18) 
ot 

^ =-iyV^P + XuVp , (19) 
ot 

for the slope u{x,t) and a canonically conjugate noise field p{x,t), replacing the 
stochastic noise ri{x,t). The field equations bear the same relation to the Fokker- 
Planck equation (|l6|) as the classical equations of motion bear to the Schrodinger 
equation in the semi-classical WKB approximation. 

To justify the weak noise limit we recall the analogy with the WKB approximation 
in quantum mechanics which, owing to its nonperturbative character, captures features 
like bound states and tunneling amplitudes, which are generally inaccessible to 
perturbation theory. Therefore, we anticipate that the present weak noise approach 
to the Burgers equation also accounts correctly, at least in a qualitative sense, for 
the stochastic properties even at larger noise strength. However, there may be an 
upper threshold value beyond which the system may enter a new stochastic or kinetic 
phase. In the one dimensional case discussed here the scaling behavior is controlled by 
a single strong coupling fixed point which can be accessed by the present weak noise 
approach. In two and higher dimension a dynamic renormalization group analysis 
predicts a kinetic phase transition at a critical noise strength (or coupling strength) 
and the weak noise approach presumably fails. 

The equations (|l^ and (|l^) derive from a principle of least action characterized 
by an action S{u' u",t) for an orbit u'{x) — > u"{x) traversed in time t 

S{u' ->u'\t) = dt dx {p^ - , (20) 

with Hamiltonian density 

n^p(^i^\/'^u + XuVu-^V^p^ . (21) 

The action is of central importance in the present approach and serves as a weight 
function for the noise-driven nonequilibrium configurations in much the same manner 
as the energy E in the Boltzmann factor exjp{—f3E) for equilibrium systems, where 
P is the inverse temperature. The dynamical action in fact replaces the energy 
in the context of the dynamics of stochastic nonequilibrium systems governed by 
a generic Langevin equation driven by Gaussian white noise. The action provides 
a methodological approach and yields access to the time dependent and stationary 
probability distributions, 

s{u' ~^u",ty 



P{u — !■ u" , t) oc exp 



(22) 



P,t{u")= lim P{u' ^u\t) , (23) 

t — >oo 

and associated moments, e.g., the stationary slope correlations 

{u{xt)u{QQ)) = iX\_du u" {x)u' {0)P{u' u" ,t)Pst{u') . (24) 
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The contact with the Fokker-Planck equation (|T6|) in the weak noise hmit is 
estabhshed by noting that (|2^) inserted in (^6|) to leading order yields the Hamilton- 
Jacobi equation dS/dt + H{p, u) — with canonical momentum p ~ 5S/Su, where the 
Hamiltonian H = J dx Ti.. 

The canonical formulation associates the conserved energy E (following from 
time translation invariance), the conserved momentum 11 (from space translation 
invariance), and the conserved area M (from the Burgers equation with conserved 
noise): 



J dxH , (25) 
U = J dx uVp , (26) 
M ^ dxu . (27) 



The field equations (|T^) and (|l^) determine orbits in a canonical up phase space 
where the dynamical issue in determining S and thus P is to find an orbit from u' to 
u" in time t, p being a slaved variable. Note that unlike dynamical system theory we 
are not considering the asymptotic properties of a given orbit. In general the orbits in 
phase space lie on the manifolds determined by the constants of motion E, 11 and M. 
Here the zero energy manifold E = plays a special role in defining the stationary 
state. For vanishing or periodic boundary conditions for the slope field the zero 
energy manifold is composed of the transient submanifolds p = and the stationary 
submanifold p = 2vu. The zero energy orbits on the p — Q manifold correspond to 
solutions of the damped noiseless Burgers equation; the orbits on the p = 2vu are 
solutions to the undamped noiseless Burgers equation with negative damping, i.e., v 
replaced by —v. In the solvable linear case of the noise driven diffusion equation for 
A = 0, i.e., the Edwards- Wilkinson equation |^^, a finite energy orbit from u' u" in 
time t migrates to the zero energy manifold in the limit t — > oo, yielding according to 
(pO|) and (|2^) the stationary distribution Pgt oc exp(— (i//A) J dxu^). This distribution 
also holds in the Burgers case and is a generic result independent of A Finally, 
in the long time limit an orbit from u' — > u" is attracted to the hyperbolic saddle 
point at the origin in phase space implying ergodic behavior in the stationary state. 
In figure we have schematically depicted possible orbits in phase space. 



In the linear case the field equations ( |18D and ( |19| ) couple p parametrically to 
u, i.e. p is slaved. In wavenumber space pk is growing, whereas uj^ driven by is a 
linear superposition of damped and growing diffusive modes, supporting the expression 
(pT|). In the nonlinear case (|l^ and (^9|) admit nonlinear soliton or smoothed shock 
wave solutions which are, in the static case, of the kink- like form in ( |l3| ) and (|l5|). 
Propagating solitons are subsequently generated by the Galilean boost (|^) and we 
recover the soliton condition (p^. 

The right hand soliton moves on the noiseless manifold p — Q and is also a 
solution of the damped (stable) noiseless Burgers equation for 77 = 0. The noise- 
induced left hand soliton is associated with the noisy manifold p = 2vu^ and is a 
solution of the undamped (unstable) noiseless Burgers equation with v replaced by 
—v. In addition the field equations also admit linear mode solutions superimposed as 
ripple modes on the solitons. The ripple modes are superpositions of both decaying and 
growing components reflecting the noiseless and noisy manifolds p = Q and p = 2i'u, 
respectively. The soliton mode induces a propagating component with velocity Xu in 
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Figure 6. Generic behavior of the orbits in up phase space. Heavy Unes indicate 
the zero energy manifold. The stationary saddle point (sp) is at the origin. The 
finite time orbit from u' to u" is attracted to the saddle point for t ^ oo. 



Figure 7. Slope field u and height profile h for the right hand and left hand 
moving kink solitons in the description of a growing interface. 



such a way that the right hand sohton acts Hke a sink and the left hand sohton as a 
source of hnear modes |^ . This mechanism wih be discussed heuristically in section 
4. In the Edwards- Wilkinson limit for A ^ the ripple modes become the usual 
diffusive modes (growing and decaying) of the driven stationary diffusion equation. In 
figure ^ we have shown the right hand and left hand solitons. 

The heuristic physical picture that emerges from our analysis now supported by a 
weak noise analysis of the Fokker-Planck equation is that of a many body formulation 
of the pattern formation of a growing interface in terms of a dilute gas of propagating 
solitons matched according to the soliton condition (|4|) with superimposed linear 
ripple modes. 



4. A growing interface 



As discussed above a growing interface can be envisaged as a many body system in the 
Landau quasi particle sense composed of matched right hand and left hand solitons 
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Figure 8. Slope field and the resulting height profile for a soliton pair 
eonfiguration. 



with superimposed linear modes. Focusing on the sohtons the right and left hand kinks 
are the fundamental growth modes corresponding to cusps in the height field; they 
are the quarks in the present formulation. They do, however, not satisfy periodic or 
vanishing boundary conditions in the slope field u; the nonvanishing boundary values 
tt+ and U- in fact correspond to a deterministic current dissipated or generated at 
the soliton centers yielding permanent profile solutions. The simplest mode satisfying 
periodic boundary conditions is the two-soliton or pair soliton configuration obtained 
by matching a right hand and a left hand soliton boosted to the velocity v = —Xu. The 
two-soliton mode has amplitude 2u and size £. By inspection it is seen that the pair 
mode is an approximate solution to the field equations (|l^) and (^9|). The correction 
terms are of the type uVu and uVp and thus correspond to local perturbations from 
a region of size ://A|u| which is small in the low viscosity limit i' 0. We assume that 
the correction can be treated within a linear stability analysis and thus gives rise to a 
linear mode propagating between the right hand and left hand solitons. This property 
is borne out by a recent numerical analysis of the field equations [^0[ . As also shown 
in the numerical analysis the pair mode forms a long-lived excitation or quasi-particle 
in the many body description of a growing interface. Subject to periodic boundary 
conditions this mode corresponds to a simple growth situation. The propagation of 
the pair mode corresponds to the propagation of a step in the height field h. At each 
revolution of the pair mode the interface grows by a uniform layer of thickness 2ui. 
In figure ^ we have depicted the pair mode and the associated height profile h. 

The soliton picture also allows us easily to understand in what sense the right 
hand soliton acts like a drain and the left hand soliton as a source with respect 
to perturbations. Considering two pair solitons superimposed on the right and left 
horizontal parts of the static solitons ( p^ and ( p^ ) it follows from (14) that for a right 
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Figure 9. Source and drain mechanism for the right hand and left hand solitons. 
The perturbation attracted and repelled by the soliton centers are modeled by pair 
solitons. 





Figure 10. Four soliton representation of the slope field and the height profile. 



hand soliton perturbations move toward the sohton center and for a left hand soliton 
perturbations move away from the soliton center. This mechanism also follows from 
the linear analysis of ripple modes . The mechanism is depicted in figure ||. 

Generally a growing interface, ignoring the superimposed linear ripple modes, can 
at a given time instant be represented by a gas of matched left hand and right hand 
solitons as depicted in figure |l^ in the four soliton case. A gas of pair solitons thus 
constitute a particular growth mode where the height profile between moving steps 
has horizontal segments. 



5. Dynamics, statistics, and scaling 



There two levels of description: the stochastic Langevin level and the deterministic 
Fokker-Planck or equations of motion level. On the Fokker-Planck level yielding the 
canonical field equations ( p^ ) and ( p^ ) the growth of the interface is interpreted in 
terms of a gas of propagating solitons (and diffusive modes) . The stochastic description 
on the Langevin level is then established in the weak noise limit A — > by computing 
the action S associated with a particular dynamical mode and subsequently deduce 
the probability distribution according to (|2^), i.e., P oc exp(— S/A). This procedure is 
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Figure 11. Dispersion law for the noise-induced left hand soliton. 



completely equivalent to the WKB limit of quantum mechanics. Here the wavefunction 
VE* and thus the probabilistic interpretation is given by ^I' oc exp{iS/h), where S is the 
action associated with the classical motion. Note that unlike quantum mechanics there 
is no phase interference in the stochastic nonequilibrium case. 

5.1. Dynamics 

The canonical phase space approach discussed in section 3 associates a formal 
dynamics with the soliton-diffusive mode gas representation of a growing interface. 
According to ( psf ) , (|2^) , and ( ^ ) the total energy, momentum and area are conserved 
in the course of the dynamical evolution of an interface on the deterministic Fokker- 
Planck level as determined by the field equations (|lj) and (|9[). Since the symplectic 
structure via the Fokker-Planck equation is associated with the noisy case only 
the noise-induced left hand soliton and, similarly, the growing linear mode, carries 
dynamical attributes. By insertion the static soliton (|l5|) thus has the energy 
E = — (16/3)Az^Uq and momentum 11 = 0. Correspondingly, the moving pair soliton 
excitation depicted in figure || is endowed with the same energy but has momentum 
n = 4i'm|u| pointing in the same direction as the velocity, it moreover conserves 
the area under propagation. Eliminating the amplitude dependence the dynamical 
characteristics of the left hand soliton is conveniently given by the soliton dispersion 
law 

1^1 -^^n3/^. (28) 

We note the fractional power arising from the nonlinear character of the soliton; the 
dispersion law is shown in figure |l^. 

5. 2. Statistics 

We consider here as a simple application the statistics of the long-lived pair soliton 
shown in figure ^. It has size £, amplitude 2m, and propagates with velocity v = —Xu. 
During a revolution in a system of size L with periodic boundary conditions the 
height field increases with a layer of thickness 2ui. Since the system is traversed in 
time t = L/v the integrated growth velocity is given by 2\v?l/L which for a single 
pair of fixed size vanishes in the thermodynamic limit. On the other hand, the local 
growth velocity dh/dt is given by 2Xu^ — (A/2)(Vft-)^ which is consistent with the 
averaged KPZ equation (|l]) in the stationary state. 
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The stochastic properties of the pair soUton growth mode is easily elucidated by 
noting that the action associated with the pair mode is given by = (4/3)j/A|M|^t. 
Denoting the center of mass of the pair mode by x we have u — v / X — x/tX and we 
obtain using ( p2| ) the transition probability 

/'(x,t)ocexp(-l^^) , (29) 

for the 'random walk' of independent pair solitons or steps in the height profile. 
Comparing (|2^) with the distribution for ordinary random walk originating from the 
Langevin equation dx/dt — rj, {'i]ri){t) — A5{t), P{x,t) oc exp(— x^/2At), we observe 
that the growth mode performs anomalous diffusion. The distribution ( p9| ) also implies 
the soliton mean square displacement, assuming pairs of the same average size, 

(x2)(i)oc(^^y^%2/. ^ (30) 

with dynamic exponent z = 3/2, identical to the dynamic exponent defining the KPZ 
universality class. This result should be contrasted with the mean square displacement 
(x^) oc At^/^, z = 2, for ordinary random walk. The growth modes thus perform 
superdiffusion. 



5. 3. Scaling 

Aspects of the scaling properties of the noisy Burgers equation are embodied in the 
scaling form ^ for the slope correlations. Here we give a set of heuristic arguments 
implying that the dynamic scaling exponent z can be inferred from the exponent in 
the soliton dispersion law (p8|); we refer to Q |l) for more details. 

Within the weak noise approximation the slope correlations are given by 
appropriate overlap integrals involving the soliton configuration |3l]. However, noting 
that the canonical weak noise formulation in general follows from a saddle point 
approximation to the Martin-Siggia-Rose functional integral, we infer that the slope 
correlation can also be expressed as the time-ordered product |32| {u{xt)u{00)) oc 
{0|T{t(xi)u(00)|0). Here the "quantum operators" u and p evolve according to the 
"quantum Hamiltonian density" ( pT] ) and |0) denotes the zero-energy stationary state. 
Displacing the field from (x,i) to (0,0), using the Hamiltonian and momentum 
operators and inserting a complete set of intermediate quasi-particle momentum states 
|n), we infer the spectral representation 

{u{xt)u{m)) (X J dU G(H) exp {Et - iUx) . (31) 

Here G(n) is an effective form factor and E and H the energy and momentum of the 
appropriate quasi particle. 

The scaling limit for large x and large t corresponds to the bottom of the quasi- 
particle spectrum and we note that only gapless excitations contribute. Assuming a 
general dispersion law with exponent (3, E oc 11^ the dynamic exponent z is given by 
the exponent f3 for the quasi-particle dispersion law. In the linear EW case the gapless 
diffusive dispersion law i5 oc H^ yields the dynamic exponent z = 2; in the Burgers- 
KPZ case the noise excites a new nonlinear gapless soliton mode with dispersion 
E oc iP^^, yielding the exponent z = 3/2. 
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6. Summary and conclusion 

We have here summarized recent work on the growth morphology and scaling behavior 
of the noisy Burgers equation in one dimension. Using a canonical weak noise approach 
to the associated Fokker-Planck equation we have discussed the growth morphology in 
terms of a gas of nonlinear solitonic growth modes with superposed linear modes. We 
have, moreover, associated the dynamic scaling exponent with the soliton dispersion 
law. 

So far the nonperturbative weak noise approach has only been implemented in the 
one dimensional case where the analysis is tractable. However, the weak noise method 
is generally applicable to generic Langevin equations driven by white noise and it 
remains to be seen whether the approach also throws light on the higher dimensional 
case. 
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